Superconductivty without inversion symmetry: MnSi versus CePt 3 Si 
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Superconductivity in materials without spatial inversion symmetry is studied. We show that in 
contrast to common believe, spin-triplet pairing is not entirely excluded in such systems. Moreover, 
paramagnetic limiting is analyzed for both spin-singlet and triplet pairing. The lack of inversion 
symmetry reduces the effect of the paramagnetic limiting for spin-singlet pairing. These results are 
applied to MnSi and CePtaSi. 
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Cooper pairing in the spin-singlet channel relies on 
the presence of time reversal symmetry (Anderson's the- 
orem); the paired electron states are related by time 
reversal and are consequently degenerated If this de- 
generacy is lifted, for example, by a magnetic field or 
magnetic impurities coupling to the electron spins, then 
superconductivity is weakened or even suppressed. For 
spin-triplet pairing, Anderson noticed that additionally 
inversion symmetry is required to obtain the necessary 
degenerate electron states^. Consequently, it became a 
widespread view that a material lacking an inversion cen- 
ter would be an unlikely candidate for spin-triplet pair- 
ing. For example, the absence of superconductivity in the 
paramagnetic phase of MnSi close to the quantum critical 
point to itinerant ferromagnetism was interpreted from 
this point of view2i4. Near this quantum critical point 
the most natural spin fluctuation mediated Cooper pair- 
ing would occur in the spin-triplet channel. However, 
MnSi has the so-called B20 structure (P2i), without in- 
version center, inhibiting spin-triplet pairing. 

Recently, superconductivity has been discovered in the 
heavy fermion compound CePtaSi, another system with- 
out inversion symmetry (P4mm)^ . The upper critical 
field H C 2 exceeds the usual paramagnetic limiting field, 
which might indicate that here nevertheless spin-triplet 
pairing is realized. Since there is no experimental infor- 
mation on the pairing symmetry in this material so far, 
it is worth examining the options for Cooper pairing in 
this case. 

The aim of this letter is to discuss two points for 
time-reversal invariant materials without inversion cen- 
ters. The first is concerned with the possible existence of 
spin-triplet pairing. The second addresses the problem 
of paramagnetic limiting (Clogston-Chandrasekar-Pauli 
limiting). The result of this discussion will be applied to 
the two materials mentioned above: MnSi and CePtsSi. 

Model: We use a single-band model with electron band 
energy £u measured from the Fermi energy where elec- 
trons with momentum k and spin s are created (anni- 
hilated) by the operators cj^ ( c ks )- The Hamiltonian 
including the pairing interaction is 

H = E Zk C ks C ks + ^EE V k,k' C ks C -ks' C -k's' C k's- 

k,s k,k' s ' s ' 

(1) 



This system possesses time reversal and inversion sym- 
metry (£u = and the pairing interaction does not 
depend on the spin and favors either even-parity (spin- 
singlet) or odd-parity (spin-triplet) pairing as required. 
Following the standard weak-coupling approach, we de- 
fine the interaction to be finite and attractive close to 
the Fermi energy with the cutoff energy e c , and to de- 
pend on the momenta only through the angular depen- 
dence. The absence of inversion symmetry is introduced 
by an additional term, H p , to the Hamiltonian which re- 
moves parity but conserves time reversal symmetry, i.e. 

-H v and TH^T^ 1 = H v . We can write such 



a single-particle term as 



' c ks c ks 



(2) 



where cr denotes the Pauli matrices and g_ k = 9 k 
(this satisfies the above condition since Icrl^ 1 = cr and 
TerT -1 = — cr). It is convenient to normalize g k so that 
the average over the Fermi surface (9 k ) k = 1, in the 
numerical calculations we will impose this constraint. We 
will keep g k arbitrary and later provide a specific form 
of g k for MnSi and CePtsSi. The normal state Green's 
function becomes, 

G°(k,iui n ) = G + (k,iu> n )a + g k ■ <xG_(fc,iw„) (3) 

where ctq is the unit matrix and 



G±(k,iuj n ) 



(iu n - e k + ) 1 ± (iw„ - e k y 



_ ( 4 ) 

e k.± =^k ±a \9k\ andfffe =9fc/ISfcl (Ifll = V^) 6 - Thc 
Fermi surface splits into two sheets with different spin 
structure. These two sheets touch whenever g(k) = 0. 

Superconducting instability: We now use the BCS de- 
coupling scheme and determine the linearized gap equa- 
tion in order to calculate the transition temperature T c : 
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A„,(fc) = -fc B T^ Yl V kik ,Gl 1 (k\ l u n )A SuS2 (k')G° s , S2 (-k\^u n ) (5) 

n,k' Sl ' S2 

The gap function is decomposed into a spin-singlet [V'(fc)] and a triplet [d(k)] part, A(fc) = {^{k)<j a + d(k) ■ a}ia y . 
For simplicity we assume that the gap functions have the same magnitude on both Fermi surface sheets. This allows 
us to write the linearized gap equations as 

V(fc) = -k B T ^fc,fc'{[G+G+ + G_G_]V(fc') + [G + G_ + G_G+]% • d(fc')} (6) 

n,k' 

and 

d(k) = -k B T ^ V k k > ^[G + G + + G-G-]d(k') + 2G-G-[g k >(g k > ■d(k')) — d(k')] + [G + G- + G-G + ]g k >'ijj(k / ) , | (7) 

n,k' 



where we have used the short notation for the products: 
G Q Gb = G a (k, iuj n )Gb{— k, — iuj n ) with a,b = ±. For 
finite a, the spin-singlet and triplet channel are coupled, 
an effect of the missing parity 7 . However, this coupling 
depends on the degree of particle-hole asymmetry or the 
difference of the density of states on the two Fermi surface 
sheets, which yields a coupling of the order a/tF *C 1. 
Thus, we ignore these coupling terms here and consider 
the "singlet" and "triplet" channel of pairing separately. 

For spin-singlet pairing we find that the transition tem- 
perature (T c ) is given by 



In 



O 



The transition temperature remains essentially un- 
changed from ksT cs = e c exp(— 1/A S ) (here T cs is T c for 
a = 0) with \ s ip(k) = N(Q)(V kk >i>(k')} k <. For triplet 
pairing the equation for T c reads 



In 



2({\d(k)\ 2 -\g k -d(k)f}f(p k )) k +0 



a 



where fcsT ct = e c exp( 
N(Q)(V k k ,d(k')) k < and p k = 

normalized gap function with (\d(k)\ 2 ) k = 1 in all 
merical calculations. The function f{p) is defined as 



(9) 

T/At) with \ t d(k) = 
a\g k \/irkBT c . We use the 



Hp) = 



n=l 



1 



1 



2n- 1 



ip 



2n- 1 



(10) 



The correction term in Eq. © suppresses T c in general. 
For a spherical Fermi surface and a — all gap functions 
with a given relative angular momentum £ have the same 
T c . Eci. ((TU|l determines how this degeneracy is lifted by 
the broken inversion symmetry. The highest T c is ob- 
tained for a state with d(k) || g k , for which the right 
hand side of Eq.@ vanishes and T c = T ct . Hence we 
conclude that spin-triplet pairing is not indiscriminately 
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FIG. 1: Transition temperature as a function of a for — 
(—k y ,k x ,0). The curves from top to bottom correspond to 



d = yk x — a 
respectively. 



yk x + xk y , and d = xk x + yk y + zk z 



suppressed in the absence of an inversion center. In prin- 
ciple, there may be spin-triplet pairing states which are 
completely unaffected by the lack of inversion symmetry, 
taking advantage of the spinor structure induced by g k . 

Structure of g: The vector ag k characterizes and 
quantifies the absence of an inversion center in a crys- 
tal lattice. In many cases the loss of an inversion center 
can be viewed as moving certain ions in the crystal lat- 
tice out of their high-symmetry position. This gives rise 
to internal electric fields that yield, through relativistic 
corrections, spin-orbit coupling 8 . Furthermore, shifted 
ions can open new hopping paths which involve atomic 
spin-orbit coupling on intermediate ions. 

We consider the form oig k for our two examples: MnSi 
and CePtsSi. We start with the space group that corre- 
sponds to the basic "point group" symmetry G. Due to 
the lack of inversion symmetry this group is reduced to 
a subgroup G'. The correction term g k ■ cr is invariant 
under all transformation of G', but not of G. MnSi has 
the cubic space group P2i. The point group is only the 
tetrahedral group T G Oh- The symmetry breaking term 
satisfying the above conditions corresponds to the irre- 
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ducible representation A 2u of Oh which maps to A\ of 
T c i. The expansion in k leads to (we assume a spherical 
Fermi surface for simplicity) 



also give a symmetry classification of the possible pairing 



9k- 



— kyk z {ky<J z 



- k z (jy) + cyclic perm. ofx,y, z (11) 



which transforms like xyz, a basis function of A 2u of Oh- 
The g-vector has 14 nodes on the Fermi surface, 6 along 
[100]- and 8 along [lll]-directions. The d-vector which 
remains unaffected by the lack of inversion symmetry is 
parallel to g k : 

d(k) = xk x (hi -k 2 z ) + yk y (k 2 z -kl)+zk z {kl-kl) (12) 

which also belongs to A 2u , has the same number of nodes 
as g and represents an /-wave spin-triplet pairing state. 

CePtaSi is a tetragonal system with space group 
P4mm. Here the removal of the inversion center leads 
to the point group C^ v £ D^h which corresponds to the 
loss of the basal plane as a mirror plane (z — > — z). One 



finds that 



9k f 



kx&v 



ky(7x 



which is a basis function 



of A 2u of D A h- This term has the form of the well-known 
Rashba spin-orbit coupling 8 . The g- vector has only two 
nodes lying along the [001]-direction. In Fig.l we show 
the reduction of various spin-triplet pairing states with 
this g-vector. The favored pairing state is of p- wave type: 
d(k) — xk y — yk x in A 2u of D^h, while other pairing 
states are severely suppressed for a > kBT c . 

With the reduced symmetry of a crystal, it is usu- 
ally not possible to find d(k) \\ g k which satisfies the 
linearized gap equation for a given pairing interaction 
Vu k ' ■ Nevertheless, we could determine a nearly opti- 
mal spin-triplet state compromising between the pairing 
interaction and the effect of H p . Very recently, Samokhin 
et al. have carried out relativistic band structure calcula- 
tions which indicate that a » ksT c in CePtsSi^. They 



states; 9 . 

The stability of the pairing state is not only decided 
by T c , but also by the condensation energy; which is de- 
termined by the shape of the quasi-particle gap in the 
weak coupling limit. The Balian-Werthamer state (with 
a nodeless gap), d(k) = xk x + yk y + zk z , is the most 
stable weak coupling state in a spherically symmetric (or 
cubic) system. However, in the presence of broken inver- 
sion symmetry, Eq.Q shows this state would generally 
have a lowered T c . Thus, for small enough a (a < ksT ct ) 
there could be a second superconducting phase transition 
below the onset of superconductivity leading to a node- 
less gap. 

Paramagnetic limiting: Lifting the degeneracy of the 
spins is detrimental to spin-singlet superconductivity, an 
effect known as paramagnetic limiting. Spin-triplet pair- 
ing is less vulnerable in this respect. In the absence of 
inversion symmetry, however, this effect of pair break- 
ing is modified. It is well-known that impurity spin-orbit 
scattering reduces the effect of paramagnetic limiting 10 . 
We show that an analogous effect occurs in systems with 
broken inversion symmetry. For simplicity, we ignore the 
effect of orbital pair breaking and include the magnetic 
field only through its coupling to the spin. Our aim is 
to demonstrate the effect of finite a on the paramagnetic 
limiting and an extended discussion for the upper criti- 
cal field H C 2 will be given elsewhere. We replace ctgy, — > 
a 9k = a 9k ~ ^ with h = hbH (note: g ^ ^ 9k)- 
The linear gap equations yield the transition tempera- 
tures for a continuous onset of superconductivity^. We 
first consider a dominant spin-singlet pairing (?/;(fc)) and 
ignore the induced spin-triplet pairing. Then the equa- 
tion determining T c is 



J 



T r: 



'k' 



a2 9\-h 2 



[{ag k + hy{ag k -hy]V* 
I 



f{p- k ) - f(j>t) 



k' 



(13) 



with p± = \ag k + h\/2irk B T c ± \ag k - h\/2Trk B T c . If 
it is possible to choose h _L g k for all k (as it is for 
CePtsSi), then in the small T c (h)/T cs limit, the para- 
magnetic limiting field obeys h' 2 \nh' = —a 2 \n(T c /T cs ); 



with h! = \h\/TrkBT cs . In particular, the paramagnetic 
limiting field diverges as T — > (Fig. 2). 

For the spin-triplet channel we obtain analogously 



hi 



\d(k)\ 2 { f{ P + k ) + f{pl) 



k' 



a 2 g% h 2 



[{ag k + hY{ag k -hYY/ 2 



f(Pp - f(Pl) 



/(pJ)-/(Pp 



\h-d(k)\ 



\g k ■ d(k)\ 



[{ag k +h) 2 {ag k -h) 2 Yn/ k 



k> 



(14) 



4 




0.0 0.2 0.4 0.6 0.8 1.0 

T<T C , 

FIG. 2: Paramagnetic limiting field for CePtsSi for different 
a (in units kBT cs ). The field is applied along the four- fold 
symmetry axis. 

For a — there is no paramagnetic limiting, provided 
d(k) ■ h = can be found for all k. According to Ea. ll4|) 
paramagnetic limiting is absent, if for all k h _L d{k) 
and d(k) \\ g^. For both the spin-singlet and spin- 
triplet cases, finite q Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) phases are often found when h ■ ^= 13 i 16 . 
The role of orbital effects on these phases is currently 
under investigation^. 

Discussion of the two examples: We start with MnSi, 
which does not show superconductivity in the vicinity 
of the quantum critical point of a ferromagnetic state. 
Given our result that spin-triplet pairing is not sup- 
pressed completely by broken inversion symmetry, it is 
useful to reexamine the reason for why spin-triplet su- 
perconductivity is not observed. As one would expect, 
the lack of inversion symmetry in this compound with 
(cubic) B20-structure is crucial. According to our analy- 
sis the pairing would have to occur in the /-wave channel 
in order to survive the spin-orbit coupling effect. The fact 
that the strongly anisotropic /-wave state is more diffi- 
cult to stabilize by a simple spin fluctuation mechanism 
than the p-wave pairing state might explains the absence 



of superconductivity in MnSi. 

Turning to CePtaSi we may adopt two different points 
of view. First, there is a protected p-wave spin-triplet 
pairing state {d(k) = xk y — yk x ). This may indeed ex- 
plain the apparent absence of paramagnetic limiting ob- 
served in polycrystaline samples^. On the other hand, 
it is important to notice that superconductivity appears 
here on the background of antiferromagnetic order (T/v ~ 
2K), and it seems more natural to assume a spin-singlet 
type of pairing. In this case, we could argue that para- 
magnetic limiting for a singlet state is rendered less ef- 
fective by the presence of spin-orbit coupling. To exam- 
ine this possibility in more detail we have determined 
the paramagnetic limiting field as a function of a using 
Eq. (|13|) for the field along the four-fold symmetry axis. 
This is shown in Fig. [3 note that this figure illustrates 
the divergent paramagnetic limiting field at low tempera- 
tures described earlier. It would be very helpful to study 
this system for single crystals, since for both the spin- 
singlet and spin-triplet cases a large anisotropy in the 
paramagnetic limiting field is predicted. The field along 
the four-fold axis should give no paramagnetic limiting 
in both cases. Moreover, the Knight shift should show 
related effects of the spin-orbit coupling. 

In conclusion, the analysis of the symmetry properties 
for the two materials MnSi and CePtsSi show that in the 
former system the effect of the lack of inversion symmetry 
leads to more severe restrictions for spin-triplet pairing 
than in the latter. Furthermore, paramagnetic limiting 
for spin-singlet superconductors is suppressed by broken 
inversion symmetry. In many respects, CePtsSi may be- 
come an ideal test system to study the effect of missing 
inversion symmetry on the superconducting phas o 6 i 7 i 12 . 
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Errata on "Superconductivty without inversion 
symmetry: MnSi versus CePt 3 Si" 

We have two corrections to make to our original paper. 
The first involves a correction in the factor an in Fig. 2. 
T he a ?fl values in Fig. 2 should be divided by a factor 

Also, we have omitted a contribution to the vector 
gfc for MnSi. Including this contribution gives = 
Q-i [k x , k y ,k z ] + 0L2 [k x (ky — k^), k y (k z — k x ), k z (k x — ky)] 
(in the above paper we have a\ =0). The vector mul- 



tiplied by a\ belongs to the A\ u representation of Oh 
and that multiplied by a 2 belongs to the A 2u representa- 
tion of Oh] both these vectors map to the representation 
A\ of T. Our result that the spin-triplet pairing vector 
dy, should be parallel to for spin-triplet superconduc- 
tivity to be stable is unchanged. Consequently, p-wave 
superconductivity is suppressed for MnSi if ct\ < a 2 - The 
relative size of ot\ and a 2 will require band structure cal- 
culations to determine. 

We are grateful S. Cunroe, A. Rosch and LA. Sergicnko 
for useful communications regarding MnSi. 



